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^ ' Abstract 



We propose to study hadron-hadron scattering using lattice QCD in an asymmetric 
box which allows one to access more non-degenerate low-momentum modes for a 
given volume. The conventional Liischer's formula applicable in a symmetric box is 
modified accordingly. To illustrate the feasibility of this approach, pion-pion elastic 
scattering phase shifts in the 1 = 2, J = channel are calculated within quenched 
approximation using improved gauge and Wilson fermion actions on anisotropic lat- 
tices in an asymmetric box. After the chiral and continuum extrapolation, we find 
that our quenched results for the scattering phase shifts in this channel are con- 
sistent with the experimental data when the three-momentum of the pion is below 
300MeV. Agreement is also found when compared with previous theoretical results 
from lattice and other means. Moreover, with the usage of asymmetric volume, 
we are able to compute the scattering phases in the low-momentum range (pion 
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three momentum less than about 350MeV in the center of mass frame) for over a 
dozen values of the pion three-momenta, much more than using the conventional 
symmetric box with comparable volume. 

Key words: 1 = 2, s-wave pion-pion scattering, asymmetric box, scattering phase 
shift. 

PACS: 12.38.Gc,11.15.Ha 



1 Introduction 



Hadron-hadron scattering experiments have offered us enormous amount of information con- 
cerning the interaction among hadrons. In these experiments, scattering cross sections and 
phase shifts are obtained experimentally in various channels with definite quantum numbers. 
On the theoretical side, although Quantum Chromodynamics (QCD) has been recognized 
as the underlying theory of strong interaction, theoretical explanation of hadronic scattering 
processes at low energies remains a challenging problem due to non-perturbative features of 
the theory thereof. Lattice QCD (LQCD) is the only systematic, non-perturbative method 
of QCD which in principle can be applied to calculate these low energy physical quantities 
from first principles using numerical Monte Carlo simulations. Calculation of hadron-hadron 
scattering phase shift is also a very important step to deepen our understanding of the strong 
interaction beyond single hadron spectrum. 

Lattice calculation of hadron scattering processes relies on a finite size method proposed by 
M. Liischer [Tll2T3ll4"ll5] in which two particle elastic scattering phase shifts (in the infinite 
volume) are directly related to the energy levels of the two particles in a finite cubic box. 
The latter can in principle be extracted in lattice simulations. Using this technique, the 
scattering length and the scattering phase shifts for pion-pion scattering in the 1 = 2, J = 
channel have been studied pi7f8|9lfT0|ll|12|13lll4|15|16j in both quenched and unquenched 
lattice QCD. There also exist lattice calculations on other hadronic scattering processes using 
various lattice actions. 



1 This work is supported in part by the National Natural Science Foundation of China (NSFC) 
under grant No. 10421503, No. 10675005, No.10575107, No.10375031, No.10675101 and supported 
by the Trans-century fund and the Key Grant Project of Chinese Ministry of Education (No. 305001) 
and KJCX3-SYW-N2 (CAS). 
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The lattice results on the pion-pion scattering phases can be compared with the experi- 
mental data and results from other theoretical methods and impressive agreements were 
seen [H|15fl6j . However, since lattice results were obtained in a finite volume, the number of 
low-momentum modes accessible for lattice simulation was limited. Part of the reason is that, 
in all previous lattice studies, hadron scattering phase shifts were calculated in a cubic box 
which has the same physical extension in all three spatial directions. In this scenario, many low 
momentum modes are degenerate in energy such as modes (1, 0, 0), (0, 1, 0) and (0, 0, 1) since 
they are related to one another by cubic symmetry. As a result, one can only access very few 
low-momentum modes in the lattice calculation in a cubic box. If one is interested in the scat- 
tering phases at more values of the scattering momenta, larger physical volumes are required 
which makes the lattice simulation very costly. In this paper, we propose to use asymmetric 
boxes to study hadron-hadron scattering on the lattice. We test this idea in a quenched study 
on the pion-pion scattering in the I — 2, J — channel. If we denote the three-momentum of 
each pion in the center of mass frame by k, within the range < k 2 < O.lGeV 2 , we are able 
to obtain scattering phase shifts at more than 12 different values of k 2 while in calculations 
with comparable cubic volumes this number is restricted to only a few. 

This paper is organized as follows. In Section 2, we briefly review the theoretical formalism for 
the computation of the phase shift in an asymmetric box, extending the finite size technique 
suggested by Luscher. The corresponding formulae are modified to the case of asymmetric vol- 
ume. Possible mixing with the J = 2 channel is discussed. In Section 3, some simulation details 
are given. Our results for the scattering phases, after chiral and continuum extrapolations, 
are then compared with known results from previous lattice calculations, chiral perturbation 
theory, dispersion relations and experimental data. Reasonable agreements are found and the 
advantage of using the asymmetric box is addressed. In last section, we will summarize this 
work and give some conclusions and outlooks. 



2 Liischer's formulae extended to an asymmetric box 

Consider a cubic box with size L x L x L and periodic boundary condition. In such a box, the 
three momentum of a single pion is quantized as: k = (27r/L)n with n = (ni, tt.2, n 3 ) e Z 3 , 
where Z represents the set of all integers. In this paper, we are interested in two-pion systems. 
Taking the center of mass reference frame of the two pions, we define k of the pion pair in a 
box as: 

= + (1) 
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where E W7T is the exact energy of a two-pion system with the two pions having three momentum 
k and — k respectively in the center of mass frame. We further define q 2 via: 

Note that, because of the interaction between the two pions in the box, the value of q 2 is in 
general not equal to n 2 with n e Z 3 . In fact, in the / = 2, J = channel, the interaction 
between the two pions is repulsive, making the value of q 2 larger than the corresponding n 2 . 

According to Liischer's method [I], two-pion s-wave elastic scattering phase shift can be 
obtained from the following formula: f 2 ] 

cot S(k) = w , 3 

where the so-called zeta function Z\ m is given by: 

n (n 2 - q 2 ) s 

In this definition, 3^m(r) = r l Yi m (Q r ), with Yi m (Q r ) being the usual spherical harmonics. 
When the physical volume is large enough, Liischer's formula ([3]) can be expanded as powers 
of 1/L. The resulting formula then relates the energy level of two hadron in a finite box to 
the hadron elastic scattering length in the infinite volume. This provides a very convenient 
way of computing scattering lengths on the lattice. The formula reads: 



Ett-jt 2772.77- - 



1 + Ci( lJ +C2 Vl7 



0(L" 6 ) (5) 



m^L 3 

where the coefficients c\ = —2.837297 and C2 = 6.375183, ao is the nn elastic scattering length. 

As explained in the introduction of this paper, many low momentum modes in a cubic box 
are degenerate in energy due to cubic symmetry. Therefore, to compute the scattering phases 
at more values of the scattering momenta, one usually has to use larger cubic volumes. This 
makes the lattice calculation more costly. In this paper, we propose to use an asymmetric 
volume which provides more non-degenerate low-momentum modes with a relatively small 
volume. Liischer's original formula ([3]) is only valid in a cubic box. To utilize similar finite 
volume techniques, we must generalize Eq. ([3]) to the case of an asymmetric box. This has 
been accomplished in [TTfTS] . 



2 This assumes that the contribution from higher angular momentum modes are negligible. In the 
cubic case, the leading contamination is from I = 4. 
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In an asymmetrical box with lattice size L x r\iL x 773L, the momentum of a single pion is 
quantized as: k = (2n/L)h with n = (77,1,77,2/772,77,3/773) and n = (m, 77,2,77,3) G Z 3 . Quantities 
k and g 2 are still defined according to Eq. ([T]) and Eq. ([5]). 

The symmetry group of the asymmetric box depends on the shape of the volume we take 
in our lattice calculation. For definiteness, we choose 772 = 1 and r/ 3 = 2 in this study and 
the corresponding basic symmetry group is D 4 which has 4 one-dimensional representations: 
Ax, A 2 , B\, £>2 and a two-dimensional irreducible representation E. Rotational symmetry is 
broken and the corresponding representations for the rotational group with definite angular 
momentum quantum numbers are decomposed accordingly: 



= Af , 1 = A 2 - + E~, 2 = Af + Bf + Bf + E~ 



(6) 



The formula for the scattering phase shifts is now modified to: 

Z 00 (l, q 2 ;r] 2 , 7/3) 



cot S(k) = m 00 (q) = 
with the modified zeta function Z\ m defined as: 



n 3/2 rj2Ti 3 q 



(7) 



n (n -<r, 

The formula for scattering length is also changed accordingly: 



0(L- 



(9) 



where the coefficients 01(772,773) and 02(772,773) can be computed once 772 and 773 are given [TT] . 
For the case 772 = 1 and 773 = 2, which is the situation studied in this paper, the two coefficients 
are found to be: 

ci(l, 2) = -1.805872 , c 2 (l, 2) = 1.664979 (10) 

Therefore, just as in the cubic case, once the two-pion energy level E^ is obtained in Monte 
Carlo simulations, the corresponding phase shift 5 can be obtained via modified Liischer's 
formula 

It is noted that the representation A\ appears in both the J = and J = 2 channel. Therefore, 
in an asymmetric box, s-wave and <i-wave scattering mix with each other. This is to be 
compared with the cubic case where the lowest mixture to s-wave is from / = 4 sector. 



3 Again, we omit higher angular momentum contributions. In this case, the leading contamination 
is from I = 2. 
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Table 1 

Simulation parameters used in this work: 






U s 


V 


Lattice 


^(Ge^" 1 ) 


Number of Confs. 


l^max 


2.080 


0.7735 


0.94 


8 2 x 16 x 40 


1.5677 


464 


0.0598 


2.215 


0.7852 


0.95 


9 2 x 18 x 48 


1.3926 


425 


0.0602 


2.492 


0.8063 


0.93 


12 2 x 24 x 64 


1.0459 


105 


0.0606 



Assuming the ci-wave scattering phases are small, one can estimate its effect on the s-wave 
phase shift as follows [T8] : 



nil 



- S (q) ~ 0(g) + o- 2 (q) tan5 2 (g) , 



fill 



where the angle 0(g) is defined via: — tan 0(g) = l/m 00 (g). The function cr 2 (g) for _D 4 sym- 
metry is given by: 

m 2 02 (q) 



°2 (q) 



(12) 



1 +moo(q) ' 

which quantifies the effect due to <i-wave mixing, On general grounds, one expects the mixing 
due to higher angular momentum to be small in the low-momentum region. To estimate its 
effect in the case of pion-pion scattering, we take experimental values for the <i-wave scattering 
phases presented in Ref. [19]. We have checked this correction to our s-wave phase shifts and 
it is found that these corrections are only of about 1 — 2%, much smaller than our typical 
error bars for the phase shifts. Therefore, in what follows, we simply neglect the effects of the 
d-wave contaminations. 



3 Simulation Details 



To test our idea of using the asymmetric box on hadron-hadron scattering, we perform a 
quenched study on the pion-pion scattering phase shift in the 1 = 2, J = channel. In this 
section, we will briefly introduce our numerical results. 



3. 1 Lattice actions and simulation parameters 



The gauge action used in this study is the tadpole improved gluonic action on anisotropic 
lattices 12011211: 



J) 

i>j 



^9 £ui ~ 36 fug " 36 ^ 



The explicit formula for the function mo2(<?) can be found in Ref. 
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^ [ 9 ul 36 u± 1 [l6) 



where Pij is the usual spatial plaquette variables and Rij is the 2x1 spatial Wilson loop 
on the lattice. The parameter u s , which we take to be the 4-th root of the average spatial 
plaquette value, incorporates the so-called tadpole improvement [22] and £ designates the 
(bare) aspect ratio of the anisotropic lattice, defined as the ratio between two spacings a s /a t . 
With the tadpole improvement, experiences show that the renormalization effects are small 
for this parameter. Thus, we have not distinguish the renormalized anisotropy and the bare 
one. The anisotropic and improvement property of the lattice action makes the calculation of 
heavier hadronic objects on coarser lattice possible. The parameter (3 is related to the bare 
gauge coupling which controls the spatial lattice spacing a s in physical units. This type of 
improved gauge action on anisotropic lattices have been extensively used in lattice calculations 
on glueballs [2Tll23T2ll25ll26] . 

The fermion action used in this calculation is the tadpole improved clover Wilson action on 
anisotropic lattice [27|28|29] whose fermion matrix reads: M. xy = 5 xy a + A xy with A given 
by: 

3 

A xy = 5 xy [l/(2K max ) + Pt^CToi^Oi + Psi^U^U + 0"23-^23 + 0"31-^3l)] 

i=l 

] (14) 

where the coefficients are given by: 

rji = v/(2u s ),rj = £/2,a= l/(2«) - l/(2K max ), 

pt = c sw (l + 0/(4^), Ps = c sw /(2ut). (15) 



In this notation, the fermion propagators with different quark masses could be solved at 
the same time usine; the so-called Multi-mass Minimal Residual (M 3 R) algorithm [301I3T1I32] . 
The bare velocity of light parameter v is tuned non-perturbatively using the single pion 
dispersion relations [291133] . The parameters K max is the largest one among all k parameters 
which corresponds to the lightest valence quark mass. The asymmetrical ratio £ is always 
fixed at £ = 5. Other parameters in our simulation are tabulated in Table [TJ 

Quenched configurations are generated using the pure gauge action (13.11) with three lattice 
sizes, 8 2 x 16 x 40, 9 2 x 18 x 48 and 12 2 x 24 x 64, corresponding to /3 = 2.080, 2.215 and 
2.492, respectively. The values of /3 are chosen such that the physical volumes for these three 
lattices remain the same. The correspondence of (3 and the spatial lattice spacing a s has been 
obtained in Ref. [33]. The physical volume of our lattices is about 5.5 fm 3 which is large 
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enough to bring the finite volume errors under control. For each set of parameters, several 
hundreds of de-correlated gauge configurations are utilized to measure physical quantities. 

3.2 Hadronic operators and the extraction of one and two pion energies 

To obtain energy levels for the single and two-pion systems on the lattice, we have to construct 
appropriate correlation functions using the corresponding hadronic operators. In this paper, 
single and two pion operators are constructed using local quark fields. For the single pion 
operators, we use: 

7r + (x, t) = — d(x, t)r 5 w(x, t), 7r~(x, t) = w(x, t)r 5 <i(x, t) 

7r°(x, t) = t)r 5 rf(x, t) - d(x, t)r 5 u(x, t)) (16) 

where w(x, t) and d(x, t) are the basic local quark field operators for the up and down quark, 
respectively. In this study, the up and down quarks are taken to be degenerate in mass so 
that isospin is a good symmetry. The operator which creates a single pion with non-zero three 
momentum k from the vacuum is obtained by Fourier transform: 

<(*) = ^£Ax,*)e-*- x (17) 

where the flavor index a of pions take values a = +, — , and V3 is the three volume of the 
lattice. By calculating correlation functions of single pion operators defined above, one can 
obtain the single pion energy at vanishing and non-vanishing momenta. 

The s-wave two-pion operators in the 1 = 2 channel are defined as: 

R 

where n labels a particular mode with three-momentum k n ; i?(k„) is the rotated three- 
momentum which is obtained from k n by applying a symmetry operation R e D 4 , an element 
of the corresponding point group. Therefore, the summation of R in Eq. ffl8l) guarantees that 
the operator thus constructed falls into the Af representation of the symmetry group whose 
continuum counterpart is s-wave for the rotational group, if the contaminations from the / > 2 
sectors are negligible. 

In order to obtain the two-pion energies, which is directly related to the scattering phases we 
want to compute, we measure the correlation matrix among different non-degenerate two-pion 
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Table 2 

The representative momentum of every mode. Mode 2, 4 and 5 have accidental degenerate modes 
which are left out in the construction of two pion operators. 



Serial Number 





1 


2 


3 


4 


5 


6 


Mode n 


(0,0,0) 


(0,0,1/2) 


(1,0,0) 


(1,0,1/2) 


(1,1,0) 


(1,1,1/2) 


(1,1,1) 


Degenerate Mode n 






(0,0,1) 




(1,0,1) 


(0,0,3/2) 





modes, using the two-pion operators defined in Eq. (fT8l) : 



C mn (t) = (Ol(t)O n (t s )) . (19) 
We then follow Liischer and Wolff's suggestion and constructed a new correlation matrix: 

n(t,t ) =c(t )-*c(t)c(t Q )-s (20) 

where to is some suitable reference time. The eigenvalue Aj(t) of this new matrix Q is: 

K(t,t ) oc e -*<«o) (21) 

It can be shown that [3] this eigenvalue avoids 0(e~ AEt ) errors and the energy eigenvalues 
could be extracted by a single exponential in t. We choose seven non-degenerate momentum 
modes to construct the correlation function matrix (fl~9l . The representative momentum of 
each non-degenerate momentum mode is tabulated in Table [2j Note that some modes might 
become accidentally degenerate in the continuum limit when rotational invariance is com- 
pletely restored. That is to say, two modes are degenerate in energy in the continuum limit 
but they are not related to one another by any D4 transformation. In Table [2] we also list these 
accidental degenerate modes. Although for finite lattice spacings, scaling violations will lift 
these degeneracies, the almost degenerate modes might make the diagonalization procedure 
unstable. Therefore, in our study of the two-pion correlation matrix, the accidental degenerate 
modes are left out in the construction of the two-pion operators. 

The single pion correlations at zero spacial momentum are constructed from the wall source 
quark propagators. Effective mass functions are then used to extract the single pion mass 
values. The mass plateaus are determined automatically by requiring the minimal of x 2 P er 
degree of freedom. All effective mass plateaus are plotted in Fig. [TJ The horizontal line seg- 
ments in these figures represent the ranges of the plateaus from which the pion masses are 
extracted. The errors for the data points are obtained from a standard jack-knife analysis. 

Similar analysis is performed for the two-pion correlation matrix. It is verified that the sym- 
metric off-diagonal matrix elements, and Cji (i 7^ j) of the matrix, are almost equal to each 
other in all cases except for some high modes at large time slices. Therefore, in our analysis, 
dj and Cji {i 7^ j) are simply averaged to construct a symmetric positive-definite matrix. In 
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Table 3 

Chiral extrapolation fits of using Scheme 1 and the continuum limit of the fitted parameter Aq 
(the corresponding value in the chiral limit). 



p 


A (GeV- 2 ) 


A l (GeV" 4 ) 


X 2 /d.o.f 


range 


2.080 


-2.30(15) 


0.44(16) 


0.05 


1-12 


2.215 


-2.27(14) 


0.59(17) 


0.04 


1-8 


2.492 


-2.15(19) 


0.82(20) 


0.09 


1-4 


continuum limit 


-2.05(36) 








X 2 /d.o.f 


0.18 









Table 4 

Chiral extrapolation of at f3 = 2.492 using Scheme 2. With others still using Scheme 1, the 
continuum limit of parameter Aq is shown. 



p 


A (GeV- 2 ) 


A l (GeV" 4 ) 


A 2 (GeV" 6 ) 


X 2 /d.o.f 


range 


2.492 


-2.24(16) 


1.15(21) 


-2.24(6) 


0.13 


1-12 


continuum limit 


-2.20(30) 










X 2 /d.o.f 


0.02 











the diagonalization procedure, we set to — 3 as the reference time. The effective energies of 
two-pion energy levels are defined as E^{i,t) = ln( ^ff^ ) where the index i represents the 
i-th eigenmode and the energy plateaus E nn are found accordingly. In Fig. El the effective 
mass plateaus at (3 = 2.492 for various modes are shown. Data for other values of (3 are 
similar. The corresponding values for k can thus be obtained from Eq. ([]]) and the values for 
S may be computed by Eq. Q for each set of bare parameters. 

3.3 Results for the scattering length 

It is found that the relative three momentum |k| for Mode is far smaller than the corre- 
sponding single-pion mass values. Therefore, SE = E W7r — 2m n can be obtained from Mode 
ignoring the small three momentum effect. It is then easy to obtain the scattering length ao 
in this channel using formula ([HD for every valence quark mass and (3 value. These results are 
then used to perform the chiral extrapolation. 

We use the quantity ao/m^, which is finite in the chiral limit, for the chiral extrapolation, as 
suggested by the CP-PACS collaboration [H]. The scale is set using the pure gauge sector with 
the Sommer scale r = 0.5fm. In Chiral perturbation Theory (ChPT), The m 2 dependence of 
the scattering length is known within Chiral perturbation theory (ChPT) [35|36f37j . However, 
it is well-known that ChPT is only effective when m n is small. The pion mass range in our 
simulation (from 0.7 GeV to 1.5 GeV) is definitely beyond the applicability range of ChPT. 
Therefore, we have attempted to parameterize our data for ao/m n using either a linear function 
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Table 5 

A summary of results for the pion-pion scattering length. 











JLQCD (LIN) [8] 


-0.0406(47) 


CP-PACS (quenched) [LI] 


-0.0558(56) 


JLQCD (EXP) [8] 


-0.0410(69) 


E865 Collaboration [H] 


-0.036(9) 


C.Liu (Scheme 1)[TOJ 


-0.0342(75) 


Current algebra [38] 


-0.046 


C.Liu (Scheme 2) [TO] 


-0.0459(91) 


CHPT (one-loop) [39] 


-0.0423(10) 


X.Du[TO] 


-0.0467(45) 


CHPT (two-loop, Roy Eq.) [37] 


-0.0444(10) 


CP-PACS (unquenched) J5] 


-0.0266(16) 


Dispersion relations |40j 


-0.0440(11) 



in 

^ = A + A im l;, (22) 

or a quadratic function in m\: 

^ = A + A x ml + A 2 m* . (23) 

These two methods will be referred to as Scheme 1 and 2. It is found that only the data at 
(3 = 2.492 show significant curvature in the pion mass regime that we are studying. For the 
other two (3 values, quadratic fits do not give statistically more favorable results. Therefore, 
we only attempted Scheme 2 for (3 = 2.492. The fitting results in both schemes are tabulated 
in Table [3] and Table HI The fittings are also illustrated in the left panels of Fig. [3] and Fig. [H 



Finally, a continuum limit extrapolation is performed to get rid of the lattice spacing errors, 
using a functional form that is linear in a 2 s . A possible linear term contamination might be 
there but the coefficient of it is too small to be visible in the fitting. In Fig. [3] and Fig. H] 
(right panels), we show the results for the continuum limit extrapolation. The straight lines 
represent the extrapolation towards the a s = limit and the final results are shown as circles. 
After the chiral and continuum extrapolations, our results for the scattering length in this 
particular channel read: 

J 0.0399(70) , Scheme 1 
a ° mn ~ { 0.0359(59) , Scheme 2 { } 

The two results are consistent with each other within errors. These results can also be com- 
pared with analogous results obtained in other theoretical calculations |38fl39|l3"7"|4T)] and the 
experiment [JT]. Our results of two schemes are both compatible with the experiment and 
the result of Scheme 1 is consistent with the results from other theoretical investigations. 
Our result for the scattering length also agrees with previous lattice results obtained by other 
groups. Table [5] summarizes all relevant results for the scattering length in this channel. 
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3.4 The scattering phase shift 



We now come to the results for scattering phases. Pion scattering phases can also be com- 
puted in the low-energy regime within chiral perturbation theory [3"51I3"(3] . However, as already 
mentioned above, the formulae thus obtained are only applicable for very light pion mass and 
low scattering momenta. We therefore used a method that has been used in quenched stud- 
ies by the CP-PACS Collaboration [T4l42|ll5] . namely we simply parameterize the scattering 
phases with a polynomial in m 2 and the momentum. Our method is a modified version of 
their methods. 

CP-PACS Collaboration [T4l42l)15j defines a scattering amplitude as follows: 

- tan£(fc) E nw 
A(m WJ k) = — (25) 

where k — |k| . Then, they used a polynomial function in both m 2 and k 2 to fit their simulation 
data. However, if some of the phase shifts data exceede the limit —90°, since the function 
A(m n , k) defined above involves the tan function which is discontinuous at —90°, this makes 
the amplitude A(m n , k) discontinuous as well which is not convenient for chiral extrapolations. 
To overcome this difficulty, we parameterize the phase shift 5 itself by a polynomial in both 
m 2 and k 2 as: 

5{ml, k 2 ) = D 00 + D 10 ml + L> 20 m£ + D m k 2 + D n mlk 2 + D 02 ¥ (26) 

The above function includes all terms with powers not larger than four. We have tried 
polynomial functions with higher powers but found that they had not improved the fitting 
quality. Note that in quenched lattice QCD, since the chiral behavior is different from true 
QCD, physical quantities can suffer from bad chiral behaviors [13]. For example, in the above 
fitting formulae, a non- vanishing constant term Dqq can exist, which would be absent in true 
QCD due to chiral symmetry. In practice, by fitting of our quenched data, we find that Dqo 
is always consistent with zero within statistical error when we regard it as a free parameter. 
Therefore, in the following discussion D 0Q is fixed to zero. 

All fitting results for Dij are tabulated in Table [6j The x 2 /d.o.f of (3 = 2.215 is somewhat 
large but still acceptable. In Fig. [61 we plot the fitted results for the phase shift S as a function 

— 9 

of momentum k while setting to zero (chiral limit). Results for three (3 values are all 
shown in this figure with different symbols. The results for different lattice spacings tend to 
agree with one another in the low-momentum limit and deviate in the large momentum limit 
as expected. We then perform the continuum limit extrapolation of various coefficients by 
using a function linear in a 2 . All extrapolations tabulated in Table [6] are good except for D 02 . 
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Table 6 

Fitted results for the scattering phase shifts 5 at each (3. The continuum limit for each parameter is 
also shown. 



p 


D 10 (GeV 2 ) 


A20 (GeV 4 ) 


An (GeV 2 ) 


D n (GeV 4 ) 


D02 (GeV 4 ) 


X 2 /d.o.f 


2.080 


-24.8(33) 


10.9(19) 


-165(11) 


-133.8(96) 


259(16) 


0.98 


2.215 


-18.8(23) 


6.7(19) 


-161.2(74) 


-104.8(61) 


190.8(99) 


2.05 


2.492 


-11.6(19) 


3.42(91) 


-191.3(92) 


-69.5(36) 


252(12) 


0.77 


continuum limit 


-1.4(41) 


-1.9(25) 


-213(19) 


-20.2(92) 


254(25) 





The continuum extrapolations are also shown in Fig. After the continuum extrapolation, 

— 2 

the results for the phase shift as a function of k are plotted in Fig. [6J with upside-down 
triangles. 



Our results can be compared with previous quenched lattice results from CP-PACS collabo- 
rations [2] and unquenched results from NPLQCD [TB]. We find that they agree with each 
other within errors. Due to the asymmetric box used in this study, we are able to compute the 
phase shifts at more values of scattering momentum compared with similar calculations using 
a symmetric box. For example, in Ref. [H], scattering phases are obtained at five values of 
k 2 in the range from 0.02GeV 2 to 0.34GeV 2 . By using an asymmetric box, even in a smaller 
range of 0.02GeV 2 to 0.12GeV 2 , we have over a dozen of data points for the phase shift which 
can be compared with results from other theoretical investigations and the experiments in 
more detail. 

Finally in Fig. [7J we have shown the the same result for the phase shifts 5 in the continuum 
limit together with the experimental results from CERN- Munich group [19]. It is seen that 
our final results agree with the experimental results within errors for k below 0.1 GeV 2 
which is about y/s = 0.6 GeV. At higher energies, our results deviate from the experimental 
results. This deviation might be caused by the systematic uncertainties in our calculation, 
e.g. quenching and chiral extrapolations. Numerically speaking, it is largely due to poor 
determination of the coefficient D 02 . However, we would like to point out that, it is clear from 
our quenched calculation that, the asymmetric volume technique advocated here would also 
be useful for unquenched studies once the unquenched configurations become available. 

4 Conclusions 

In this paper, we propose to study hadron-hadron scattering processes on lattices with asym- 
metric volume. This setup has the advantage that it provides much more non-degenerate 
low-momentum modes with a relatively small volume, allowing more detailed comparison 
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both with the experiments and with other theoretical results. To illustrate the feasibility of 
this proposal, pion pion scattering length and scattering phases in the / = 2, J = channel are 
computed within quenched lattice QCD using clover improved lattice actions on anisotropic 
lattices. Our quenched results indicate that the usage of asymmetric volumes indeed allow us 
to access much more low-momentum modes than in the case of cubic volume of similar size. 
For k 2 in the range of 0.02GeV 2 to 0.12GeV 2 , we have over a dozen of data points for the phase 
shift, much more than that in the cubic case with similar volume. It is also noted that, in 
the low-momentum region, after the chiral and continuum extrapolations, our results for the 
scattering length and the scattering phase shifts are in good agreement with the experimental 
data and are consistent with results obtained using other theoretical means. 

Although our calculation is now performed in the quenched approximation, similar calcula- 
tions are also possible in the unquenched case once the gauge field configurations become 
available. Finally, we have only computed scattering length and phases in the 1 = 2, J = 
channel. Phenomenologically speaking, other channels, in particular I = J = channel, are 
more interesting. However, this channel is difficult for two main reasons: One needs to per- 
form a full QCD calculation otherwise the theory is sick in the chiral limit [13]; one has to 
deal with vacuum (disconnected) diagrams which significantly increase the amount of com- 
putational cost. Also interesting and equally challenging is the I = J = 1 channel where one 
would expect to see a rho resonance [H]. We expect the use of asymmetric volumes should 
also be beneficial in these studies since a lot more low-momentum modes become accessible 
in an asymmetric box. 
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p=2.215 




Fig. 1. The single pion effective mass plateaus are shown. The horizontal line segments in the figure 
represent the fitting ranges of the plateaus. 
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Fig. 2. Effective mass plateaus of two-pion system at (3 = 2.492 after the diagonalization procedure. 
The horizontal line segments in the figure represent the fitting ranges of the plateaus. Results for 
other values of (3 are similar. 
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Fig. 3. Left panel: chiral extrapolation of the quantity using Scheme 1 as described in the paper. 
Right panel: the corresponding continuum limit extrapolation of the quantity The extrapolated 
result is indicated by a solid circle near a s = 0. 




m K 2 (GeV 2 ) a s 2 (fm 2 ) 



Fig. 4. Left panel: chiral extrapolation of at (3 = 2.492 using Scheme 2 while keeping other beta 
values extrapolated in Scheme 1. Right panel: the corresponding continuum limit extrapolation of 
the quantity ■^P-. The extrapolated result is indicated by a solid circle near a s = 0. 
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Fig. 5. The continuum limit extrapolation of the parameter Dij. The solid circles near a s = 
represent the corresponding continuum limit values. 
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Fig. 6. The values of the phase shifts 5 after the chiral extrapolations at three values of (3. The data 
points are labeled using squares, circles and triangles for three values of (3 respectively. Also shown 
are the continuum limit results of the phase shift which are labeled by upside-down triangles. 




Fig. 7. Comparison of our lattice results for the scattering phase shifts with the experimental data 
from CERN- Munich [19]. Results are consistent with each other for k 2 below 0.1 GeV 2 which roughly 
corresponds to the center of mass energy of about 0.6GeV. 
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